Wind instability of a foam layer sandwiched between the atmosphere and the ocean 
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The Kelvin-Helmholtz (KH) instability of short gravity waves is examined in order to explain the 
recent findings of the decrease in momentum transfer from hurricane winds to sea waves. A foam 
layer between the atmosphere and the ocean is suggested to provide significant stabilization of the 
sea-water surface by the wavelength shift of the instability towards smaller scales. It is conjectured 
that such stabilization leads to the observed drag reduction. The problem of a three-fluid system 
with large differences in densities provides an extension to the fundamental KH problem in fluid 
mechanics. 
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Introduction. — Results of direct measurements extrap- 
olated from weak to strong winds predict a linear increase 
of momentum transfer from wind to sea waves. The 
present study is motivated by recent findings of satura- 
tion and even decrease in the drag coefficient (capping) 
in hurricane conditions that is accompanied by produc- 
tion of a foam layer on the ocean surface A possi- 
ble explanation for the phenomenon is the development 
of a foam layer at the air-sea interface. The principal 
role of such an air-water foam layer in energy dissipa- 
tion and momentum transfer from hurricane wind to sea 
waves has been first suggested in Q. Winds generate 
waves on the ocean surface with a wide spectrum of wave 
lengths. The longest waves, hundreds meters of length, 
attempt to catch up with the wind, while the steeper 
short waves break out and play a dominant role in drag 
production (ij. When the wind speed exceeds storm 
force ( 24m/s), wave breaking creates streaks of bubbles 
near the ocean surface. As the wind exceeds the hurri- 
cane force ( 32m/s), streaks of bubbles combined with 
patches of foam cover the ocean surface. When the wind 
speed exceeds 50m/s, a foam layer completely covers the 
ocean surface 

Nowadays, there is a little hope for a comprehensive 
numerical calculations of the drag coefficient reduction 
that includes a detailed description of the wave break- 
ing and foam layer production. Indeed, up to now there 
is no complete understanding of the phenomenon. In 
the present study, the intermediately short wave Kelvin- 
Helmholtz instability (KHI) of a foam layer be- 
tween the atmosphere and the ocean is investigated in 
order to qualitatively explain the drag reduction phe- 
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nomenon. Such three-layer system exhibits a high con- 
trast in densities of constituting fluids p a -C Pf *C pw 
The present study is not concerned with the formation 
mechanism of the foam layer by the hurricane but rather 
focuses on how a foam layer isolates the lower atmosphere 
from the sea surface The existence of the foam layer on 
the ocean surface is postulated and supported by obser- 
vations (see [lj], Q and references therein). The present 
modeling demonstrates a new effective mechanism to sta- 
bilize the sea surface by a thin foam layer between the 
atmosphere and the ocean. However, beyond that par- 
ticular application, the current work addresses a funda- 
mental problem in fluid mechanics which provides a gen- 
eralization of the classic KHI. Thus, the peculiarities of 
three-layer systems with large differences in the densities 
may be of interest to a wide range of applications in the 
laboratory as well as in geophysics and astrophysics. 

The physical model. — A piecewise constant approxi- 
mation for the equilibrium densities and for the longi- 
tudinal velocities of the water, foam and air pj and Uj, 
(j = a, /, w) is employed: 

p = p w , U = U l0 = for y < 0, 
p = Pf,U = Uf for < y < Lf, 

p = p a ,U = U a for y > 0. (1) 

Here U a is the known constant velocity of the wind, while 
the constant foam layer thickness Lf and velocity Uf 
are the widely unknown parameters of the foam layer in 
hurricane conditions. In addition, it is assumed that the 
equilibrium state is in hydrostatic equilibrium, namely, 
dPj/dy = —gpj (g is the gravity acceleration). 

The equations of motion that govern the dynamics of 
the system in each of the three layers, and the appro- 
priate boundary conditions are applied at the foam layer 
interfaces. The equilibrium state is perturbed as follows: 

$(x,y,z) = F(y) + F'(x,y,t), (2) 

where $ stands for any of the physical variables, and F 



and F 1 denote the equilibrium and perturbed values, re- 
spectively. The latter are assumed to be of the form 
F' = f'{y)^xp{—iujt + ikx) with real k and complex 
u> = u> r + i^i- Thus, the amplitudes /' that satisfy the 
boundary conditions at y = ±00 are given by: 

fa = fa exp(-fcy), f' w = f w exp(fcy), 

// = /-/ exp(-fcy) + / +/ exp(ky), (3) 

where tilde denotes constant magnitudes. 

Finally, capillary and viscosity effects are neglected for 
both the equilibrium and perturbed states (see the sec- 
tion Results and discussion) . Substitution of Eqs. (l)-(3) 
into the linearized Euler equations and applying the con- 
tinuity conditions of normal velocity and pressure at the 
foam interfaces, yields the quartic dispersion relation for 
phase velocity C Q : 

2(H a + H w ) + (E- l)(H a + 1)(H W + 1) = 0, (4) 

where 

p w (U w - C) 2 - (p w - pf)g/k 



C — — , H u 
k 



E = exp(2kL f ),H a 



PfiUf-Cf 

Pa{Ua - C? - (pf - Pa )g/k 

p } {Uf-cy 



(5) 



Before turning to the study of the foam layer effect, it is 
noticed that in the limit Lf = 0, or, cquivalently. cither 
Pf = p w or pf = p a . Eq. (4) is reduced to the classic 
dispersion relation H a + H w = for KHI [B| : 

pw(k U w - lu ) 2 + p a (k U a - w ) 2 = k a g(p w - p a ), (6) 

where the subscript denotes the foam-free parameters. 

Asymptotic analysis. — First, the limit of low air-water 
density ratio, p a /p w — e 2 *C 1, (e 2 ~ 10" 3 ) is applied 
to the classic two-layer case described by Eq. (6) with 
U w = 0, in order to obtain an estimate for the various 
physical parameters: 

lu = yjgko - e 2 k 2 U 2 L + 0(e 2 k Q U a , egk Q /uj Q ). (7) 

Doing so, it can be concluded that the classic two-fluid 
KHI is excited in the short wavelength regime: 

fc £* - k*L* = 1/e 2 , w £*/t/* ~ 1/c, C /U* ~ e, (8) 

where U* = U a , L* = U 2 /g, while the superscript as- 
terisk denotes the marginal values of the parameters. 

Back to the general case of three-fluid systems, it is 
assumed that the water content in the foam, a w , is small 
(low water content is a characteristic feature of air- water 
foams). As a result, a w ~ 0.05 is scaled with e and yields 



Here = p w , pf = a a p a + a w p w , a a = 1 - a w . Assum- 
ing now that the three-fluid system operates in the same 
regime that gives rise to the KHI in the classic air-water 
system, the following scales are adopted: 



1 luL* 



c_ 



(10) 



Further assuming that the foam layer thickness is much 
less than the characteristic length, Lf/L* <C 1, (L* ~ 
250m for U a ~ 50m/ s), while the foam velocity is much 
less then the wind velocity and much larger the phase 
velocity e - C/U* < Uf/U» < 1: 

Uf/U*~e a ,L f /L«~e\ < a < 1, < 6, (11) 

which yields the following estimates for Eq. (4): 

H a ~ H w ~ e 1 " 2 ", E ~ exp(e b ~ 2 ). (12) 



Inserting the scaling (12) into Eq. (4), and applying 
the principle of the least degeneracy Q of the three-fluid 
problem, results in a = 1/2, 6 = 2, which means: 



U f ^ £ l/2 L f . ^0 



Pa 
P* 



(13) 



where Xq = 27r/fcg. Following relations (13), the wave 
number and frequency are rescaled as follows: 

k = k/k*~e°, lo = u/^/g~k* - e°. (14) 

This yields the dispersion relation to leading order in e: 



i 



2(k- k 2 ) -IE- l){k-'K, - k){KJ l + 1) 

1 , (15 

2 + {E-l){Kj 1 + 1) V ' 



where E = exp(2kLf), while the rescaled foam thickness 
Lf, and the equilibrium ratio of the foam-to- air dynamic 
pressure Kf (0 < Kf < 1) are given by: 



Lf = k* L f ^e°, K f = 



PaU 2 



(16) 



Thus, the system stability is parameterized by the dimcn- 
sionless foam velocity and thickness or, equivalently, Kf 
and k^Lf, which has a meaning of a bulk foam Richard- 
son number Rif scaled by p a /pf = e 2 /oc w ~ e: 



Rif = k* Lf, Rif = —g 



Pf AC/2 pf 



where AU = U a — U w = U* and Ap = p a — p w « — p*. 

Two particular limits of Eq. (15) are readily obtained, 
namely, the foam- free limit (H w + H a = for Lf = 0) : 
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and the foam-saturated limit (H w + 1 = for Lf = oo): 

vk =i ^-^ hf=KLf=o °^ (18) 

which differs from Eq. (17) by replacing fcg, c^o with 
k*oo = K/Kf> (0 < Kf < 1). Comparison of 
these two limits demonstrates the stabilizing effect of 
the foam due to the decrease of the marginal wavelength 
from the foam free Aq = 27r/fcg to the foam-saturated 
A^, = 2Tr/k^ XJ value. The growth rate LUi decreases from 
the foam-free u>io to the foam-saturated w^oo value. The 
definition for fc^ = k^/Kf is used in order to express 
K f through A^: Kf = fc 7*&> = A^/^Tre 2 ^). The 
intermediate wavelength value A « A^ « lm is chosen 
for further estimations from the wavelength range of the 
drag responsible waves ~ 0.1 — 10m. In turn, a typical 
height h « 0.1m is expressed from Stokes heuristic rule 
for the critical steepness of breaking waves [l(| ■ Con- 
sequently, the value Kf ss 0.5 is adopted that results in 
Uf = tU a \J K f I a w ~ 5m/s at a w w 0.05, U a ~ 50m/s. 
Figure Q] depicts the growth rate as a function of the 



C0 ; 




FIG. 1: Growth rate u>i = cji/^/gkg vs wave number, k = 
k/kg, for the typical foam-layer thicknesses, Rij = Lf = kgLf 
and the ratio of the foam/air dynamic pressure Kf = 0.5. 

wavenumber. As can be seen, the growth rate decreases 
as the foam layer thickness is increased and approaches 
its saturated limit already at k^Lf ss 1. The depen- 
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FIG. 2: Growth rate Wj = tjJ%/ -J gkg vs foam-layer thickness, 
Rif=Lf = kgLf, for the typical wave number k — k/kg and 
the ratio of the foam/air dynamic pressure Kf = 0.5. 

dence of the growth rate u>i / y^gk^ on the foam-layer 



thicknesses is depicted in Fig. [2] For sufficiently short 
waves (fc/fcp > l/Kf) the growth rate strongly drops 
from the foam free value at k^Lf=0 to its saturation level 
at foam-layer thickness k^Lf sa 1. The growth rates of 
perturbations with longer waves (fc/fcg < l/Kf) sharply 
decrease with the increase k^Lf, till total stabilization 
at a finite value of k^Lf is achieved. These two cases 
are separated by the threshold curve (fc/fcg = 1/Kf) f° r 
which the growth rate vanishes at k^Lf >> 1. 

The marginal wave number k* satisfies the eigenvalue 
equation for the three-layer system: 

exp(2k*L f ) = 1 2 .yfif - (19) 

V /; l + KJ 1 l-k*/^ v ' 

As in the classic two-fluid system, to leading order in e, 
the waves propagate with phase velocity C = ui/k with- 
out amplification for k/k* < 1, and amplify with zero 
phase velocity for k/k* > 1. The value k* monotonically 
increases with k^Lf from the foam-free value k* = fcg to 
the foam-saturated value k* = k^ = k^/Kf. 

Results and discussion. — The atmosphere-ocean in- 
teraction in hurricane conditions creates a foam layer be- 
tween the atmosphere and the ocean. This provides for 
an effective mechanism of the sea surface stabilization. 

The analysis of the KHI is treated asymptotically in 
two small parameters: air-water density ratio ~ e 2 and 
water content in the foam ~ e. The system stability 
is parameterized by the dimensionless foam velocity Uf 
and thickness Lf (or, equivalcntly, the dynamic pressure 
ratio Kf and Richardson number Rif). Due to lack of 
observations or modelling data in hurricane environment, 
they are first estimated as -L//L* = e 2 and Uf/U* ~ s/e 
by applying the asymptotic principle of least degeneracy 
of the problem. Then w 0.25m at [/» w 50m/ s 

is evaluated by the condition that the growth rate ap- 
proaches its minimal saturated value at Lf = e 2 L*, 
and further increase Lf is ineffective, as if the foam layer 
is of infinite thickness. The value Lj' is of the or- 
der of the experimentally observed values (0] and refer- 
ences therein). The single fitting parameter of the model 
Kf = A^./Aq ~ 0.5 (Uf ~ 5m/ s) had been estimated 
through an intermediate value of length of drag respon- 
sible waves (A w lm). The value of the wavelength ratio 
exhibits the instability shift towards smaller wavelength 
scales. Thus, the foam layer reduces the foam-free wave- 
length Aq approximately by a factor 2 to the foam satu- 
rated limit A£o already at Lf m L^^ . This scale-down in 
the characteristic unstable length scales provides a qual- 
itative link between the linear stability modeling and the 
role of the foam layer in the air-sea momentum exchange. 
To see that, the local correlation, based on the dimen- 
sional grounds, z/X = F(h/X) between the ocean surface 
roughness z/X and the wave steepness h/X, is examined 
in a vicinity of intermediate values of height and length of 
drag responsible breaking waves. It is similar to the cor- 
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relation ll| for pre-hurricane conditions, but with the 
local values h, A instead of significant wave height and 
peak wavelength. Noting that the breaking process does 
not completely destroy the waves, but rather tears off 
their tops, when their steepness exceeds a critical value ( 
determined by nonlinear effects), the complete covering 
of the ocean surface by the foam layer occurs when the 
critical steepness f/fO [T3| is achieved for drag responsi- 
ble breaking waves. As a result, the roughness is reduced 
along with the wave length by a factor Aq/A^ = KJ 1 ps 2 
(for A^ ps A ps lm, h ps 0.1m, Kf ps 0.5) due to the foam 
effect. Remarkably, such simple scalings are supported 
by the observations of the roughness and drag reduction 
presented in [lj. The results are physically transparent, 
since in the foam saturated system the foam layer totally 
separates the air flow from the sea surface, and the three- 
fluid system becomes close to a two-fluid foam- water sys- 
tem. Formally this corresponds to substituting the foam 
density and velocity instead of those parameters for the 
air in the classic two-fluid model. 

Finally, the main assumptions adopted in the present 
study are discussed. First, it is noted that the assump- 
tion of foam thickness uniformity, breaks down due to 
the foam accumulation in troughs of high and long ocean 
waves (seen in a photograph of the sea surface before its 
complete coverage by a foam [lj). The drag of the hurri- 
cane induced surface waves comes mainly from intermedi- 
ately short waves (of ~ 0.1 — 10m [ll])- Thus, although 
the breaking waves under consideration (of ~ lm) are 
strongly modulated by long waves (of ~ 10 2 m), such a 
slow variation of the foam thickness may be taken into ac- 
count in the next approximation. It is assumed that the 
air-sea exchange of the highly solvable carbon dioxide or 
oxygen leaves the equilibrium water content to be small 
in compliance with the determining feature of gas-liquid 
foams. Zero compressibility and viscosity approximation 
is commonly accepted in KHI of air-water systems 
[f| . The foam compressibility may be ignored within the 
same accuracy as the air one. Indeed, using the smallness 
of the air Mach number M a = U a /C a and noting that the 
foam-to-air sound velocity ratio Cf/C a ~ ypa/Pf ~ v 7 ^ 
[l3| is of the same order as Uf/U a ~ y/e (see Eqs. (13)), 
it is obtained that M f = u//C f ~ M a < 1. Although 
the foam viscosity data in hurricane environment is un- 
available, artificial foam viscosities are known to be sig- 
nificantly larger than the viscosity of its liquid and gas 
constituents. On the other hand, natural sea foams are 
expected to have lower viscosity than their artificial coun- 
terpart due to lack of man-made surfactants and a larger 
effective size of the foam bubbles (of~ 0.2 — 2mm jij). In 
any case, the stability behavior regarding the growth rate 
of the shear viscosity (see e.g. [H|]) ignored in the present 
study can only enhance the foam stabilizing effects in the 
range of the intermediately short waves. Ignoring the 
capillary effects is valid at the water-foam interface since 
the foam is composed of the same sea- water. At the air- 



foam interface the value of the surface-tension coefficient 
may be naturally assumed equal to the surface tension 
between air and sea-water. Since the latter is much less 
than the surface tension between air and fresh-water due 
to the effect of the surfactants, its influence on KHI is 
rather small for the short waves (of ~ lm). The bubbly 
liquid, spray and foam coexist in hurricane environment. 
The high contrast in three-fluid densities is the principal 
feature of the foam-layer system p a <C p/ <C p w , dis- 
tinguishing it from the layers of bubbly liquid or spray, 
with p\y ps p w or p s ~ p a , when the system will be close to 
the two-fluid air-water configuration. Indeed, the bubbly 
liquid superposed the water hard to be distinguished due 
to a small contrast in densities (p ps p w , since a w ~ 1), 
and the bubbly liquid layer may be dropped from the KHI 
study. The spray between the air and the foam is well 
distinguished from them by a contrast in densities. In- 
deed, typically pj ss a w p w , and p s = a w p w + a a p a ps 2p a 
at a w ps 0.001 for spray and a w ps 0.05 for foam, and 
p s ps 2p a « Thus, the spray may be described by 
a more detailed velocity and density air profiles within 
the present air-foam-water configuration. This however 
will increase the system uncertainty, since the spray layer 
thickness and velocity arc rather unknown. 

Helpful discussions with A. Soloviev and V. 
Chernyavski are gratefully acknowledged. 
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